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We present large scale electronic structure based molecular dynamics simulations of liquid methane 
at planetary conditions. In particular, we address the controversy of whether or not the interior of 
Uranus and Neptune consists of diamond. In our simulations we find no evidence for the formation 
of diamond, but rather sp 2 -bonded polymeric carbon. Furthermore, we predict that at high tem- 
perature hydrogen may exist in its monoatomic and metallic state. The implications of our finding 
for the planetary models of Uranus and Neptune are in detail discussed. 

PACS numbers: 31.15.aq, 71.15.Pd, 96.15.Nd, 96.15.Pf 



Being the most abundant organic molecule in the uni- 
verse, liquid CH4 at high temperature and pressure is 
of great relevance for planetary science. The here con- 
sidered pressure and temperature conditions follow the 
isentrope in the middle ice layers of Neptune and Uranus 
at a depth of one-third the planetary radius below the 
atmosphere. The gravity fields and mean densities of the 
outer gas giants Neptune and Uranus allude to a three- 
layer model: a relatively small central rocky core com- 
posed of iron, oxygen, magnesium and silicon, followed 
by an ice mantle and a predominantly hydrogen atmo- 
sphere. The middle ice layer consists of CH4, NH3 as 
well as H 2 and, in spite of its name, is not solid but 
gaseous in the outer atmospheres and a hot liquid in the 
interior. At variance to the planetary models of Saturn 
and Jupiter, the observed values for mass and radius in- 
dicate that hydrogen cannot be an integral part of either 
Neptune and Uranus. Since it is moreover not primordial, 
the detected abundance of hydrogen in the atmospheres 
of both planets implies that it may initially originate from 
deep within the planets and brought to the outmost layer 
by convection, where it does not substantially contribute 
to the total mass 

In any case, information on the interior structure of 
Neptune and Uranus are scarce and experimentally only 
indirectly accessible by means of Voyager II flyby mea- 
surements [H, shock- wave compression 0, SLas well 
as laser- heated diamond anvil cell experiments [6|. Even 
though CH4 is the most stable hydrocarbon at ambi- 
ent conditions, based on these shock-wave experiments 
as well as theoretical ground state calculations jjj], it 
has been suggested that CH4 may dissociate around 
P = 20 GPa and T = 2000 K into H 2 and diamond @|. 
While there is little doubt that in the cores of Uranus 
and Neptune CH4 dissociates into diamond, this would 
be anyhow rather consequential as it implies that in the 
interiors of these giant planets there is no CH4 at all, but 
a huge diamond mine instead. 



On the other hand, ab-initio molecular dynamics 
(AIMD) simulations predicted that the formation of di- 
amond is preempted by the appearance of hydrocarbons 
Notwithstanding that their finding had been sub- 
sequently confirmed by laser-heated diamond anvil cell 
experiments, which, at pressure P = 19 GPa and tem- 
perature T = 2000 — 3000 K, indicate the presence of 
both polymeric carbon as well as diamond @ . This view 
was further strengthened by subsequent AIMD simula- 
tions, even though none of them found any evidence for 
diamond formation 
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Nevertheless, AIMD simula- 
tions are particularly appropriate to directly probe CH4 
under the extreme pressure and temperature conditions 
predominating in the middle ice layer, in particular as 
here in either case covalent bonds are broken and formed. 
Moreover, all of the AIMD simulations show that the 
intricate interplay between temperature and pressure is 
essential to grasp CH4 at planetary conditions, where co- 
valent C-H bonds are broken by heat, while compression 
favors condensation of the dissociated carbon atoms. It 
is therefore suggestive that in AIMD simulations at even 
higher pressure, but still in the middle ice layer, carbon 
may nonetheless spontaneously transform into diamond. 
However, what causes the large discrepancy in the pres- 
sure between theory and experiment when diamond is 
formed is unknown. 

In this Letter we revisit the behavior of liquid CH4 
by means of a novel field-theoretic approach to linear 
scaling AIMD simulations However, contrary to 

previous AIMD simulations, where the considered sys- 
tems sizes have been rather small 

H H E3, here we 
use as many as 1000 CH4 molecules in a periodic cu- 
bic simulation box of length L = 25.55 A as our unit 
cell. All of our calculations have been performed in the 
canonical ensemble at T = 2000 — 6000 K and volume 
V = 10.04 cm 3 /mol, which corresponds to the second- 
shock at P = 92 GPa and T = 4000 K of a two-stage 
light-gun shock compression experiment Q ■ 




Figure 1. (Color online) Comparison of the C-C PCF at 2000 Fi g ure 2 - ( Color online ) The C - H PCF at 2000 K ( solid red 
K (solid red line), 4000 K (dashed blue line) and 6000 K (dot- line )> 4000 K (dashed blue line) and 6000 K (dot-dashed black 
dashed black line). line )- 



For all of our large scale MD simulations we have em- 
ployed the self-consistent tight-binding model for hydro- 
carbons as implemented in the CMPTool code fl7l | . 
The atoms are propagated by a modified Langevin equa- 
tion using a discretized time step of At = 0.5 fs 1^, 
Since we are dealing with a large disordered system at 
finite temperature, the Brillouin zone is sampled at the 
r-point only. 

Well-equilibrated and long trajectories are essential to 
ensure an accurate sampling. To that extend we have at 
first carefully equilibrated each of our our simulations at 
T = 2000,4000, and 6000 K, before accumulating statis- 
tics for overall 50 ps. Even though dissociation processes 
typically happen on rather short timescales, it is impor- 
tant to note that the temperature for dissociation and de- 
hydrogenation as determined by direct MD simulations 
only represents an upper bound. A major advantage of 
our novel grand-canonical simulation technique is that, 
at variance to conventional ground-state AIMD simula- 
tions, excited electrons can be employed. Due to the 
fact that they are known to dramatically weaken covalent 
bonds J20j, and therefore may facilitate the dissociation 
of CH4, we have hence chosen to be identical with 

the nuclear temperature. Nuclear quantum effects, such 
as zero-point energies, are less important for the high 
temperature regime examined here and are therefore ne- 
glected. However, entropy effects have been shown to be 
very relevant, so that the dissociation of CH4 is suppos- 
edly much more sensitive to temperature than it is to 
pressure 
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The partial pair correlation functions (PCF) of our 
simulations are shown in Figs. QJ31 As can be seen in 
Fig. [U as well as Fig. El at T = 2000 K essentially no 
covalent C-C and H-H bonds are present. The only re- 
maining significant peaks at 2.81 A and 1.74 A represent 



the average C-C and H-H distances between two adjacent 
CH4 molecules, respectively. The insignificant peaks at 
1.47 A in Fig. [1] and 0.75 A in Fig. [3] do not point 
to an onset of dissociation, as proposed by experiment 
[(J || , but are rather due to fleetingly broken C-H bonds 
caused by finite temperature. Consequently, the sharp 
intramolecular peak in Fig. [5] at around 1.075 A can be 
ascribed to covalent C-H bonds. The corresponding coor- 
dination numbers, as obtained by integrating the associ- 
ated PCF's up to their first minima, are shown in Tablefl] 
In the case of C-H the partial coordination number is 
3.984, which indicates that the liquid at T = 2000 K, ex- 
cept for single fleetingly broken C-H bonds, is nearly ex- 
clusively made up of undissociated CH4 molecules. This 
view is consistent with other AIMD studies 0, [ill 12 1 , 
but at variance to theoretical ground-state calculations 
01, as well as experimental measurements 0, H no 
signature for dissociation have been found. 

For the most relevant case at T = 4000 K and P « 
100 GPa the situation is much different and evidences 
for dissociation can indeed be observed. As can be seen 
in Fig. [1] covalent C-C bonds are appearing just like co- 
valently bonded H2 dimers, as shown in Fig. [3J As a 
consequence, the height of the intramolecular C-H peak 
in Fig. [2] is much reduced, though still existing. From 
Table U it can be deduced that nearly half of the covalent 
C-H bonds are broken, which indicates that methane does 
dissociate only partially to form hydrocarbon chains with 
mainly two and three carbon atoms, as well as H2. More 
precisely, the CH4 molecules dissociate and recombine to 
form C2H6 and to a smaller extend CaHs, which is in 
agreement with previous AIMD studies |9j, |l0|, |l2[ . How- 



ever, we find no sustained signature for the presence of 
C2H2 , which has been detected in the atmosphere of Nep- 
tune |2l|. But, we do find seeds of somewhat longer sp 2 - 
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Figure 3. (Color online) The H-H PCF at 2000 K (solid red 
line), 4000 K (dashed blue line) and 6000 K (dot-dashed black 
line). 



bonded chains and ring-like carbon structures, but defi- 
nitely no signature of sp 3 carbon bonds, i.e. no diamond- 
like carbon. This is consistent with the computed vibra- 
tional density of states of Spanu et al. (HI, who report a 
noticeable feature at 1600 cm -1 that can be attributed 
to threefold coordinated carbon atoms in graphite-like 
configurations. That is to say, that on the one hand 
our calculations are in agreement with experiment by im- 
plying that in the middle ice layer CH4 molecules itself 
are not present, but merely its dissociated constituents, 
which confirms that the interior chemistry of Uranus and 
Neptune is more complex than previously assumed. On 
the other hand, our results differ in the sense that at 
T = 4000 K we do not find any evidence for diamond- 
like carbon, as reported by the very same experiments 

SSI. 

Even deeper within the planet at even higher tempera- 
ture of T = 6000 K, the remaining CH4 have fully disso- 
ciated as indicated by the vanishing intramolecular C-H 
peak in Fig. O On the other hand, the first peak in 
Fig.[TJ which is due to C-C bonds, as well as the covalent 
H-H peak in Fig. are even more pronounced as is the 
case for T = 4000 K. As shown in Table Q] the partial 
C-H coordination number is rather small, which entails 
that contrary to T = 4000 K hydrocarbon chains are no 
longer present, but have completely dehydrogenated into 
polymeric carbon and hydrogen. As before, no evidence 
for sp 3 -bonded diamond could be found, which indicates 
that even higher pressures are required to condense car- 
bon into diamond. 

However, in contrast to previous AIMD simulations 
[I ES G3 we find for the first time that at T = 6000 K 
hydrogen is no longer solely molecular, but a noticeable 
fraction of monoatomic hydrogen can be identified. This 
immediately suggest that at T = 6000 K the present 



Table I. Partial coordination numbers, as obtained by inte- 
grating the associated PCFs up to their first minima, for all 
investigated temperatures. 



Temperature 


C-C 


C-H 


H-H 


2000 K 
4000 K 
6000 K 


0.080 
1.244 
2.556 


3.984 
2.298 
0.417 


0.004 
0.408 
1.115 



hydrogen molecules are on the verge of a liquid-liquid 
phase transition into an atomic fluid ph ase, which is 
in agreement with recent calculations 22j |23j- Due to 
the fact that dissociation is believed to be accompa- 
nied with the metallization [24l 25 1, our prediction of 



liquid atomic hydrogen may lead to an explanation for 
the large magnetic fields of planets such as Uranus and 
Neptune through a dynamo-like mechanism by electrical 
currents in the liquid metallic regions of their interiors 
[26| . Since we find that liquid CH4, where it is stable at 
T = 2000 K, is a wide band-gap insulator, it can only 
fully dehydrogenated contribute to the magnetic field in 
the form of metallic hydrogen, which indeed has been 
established experimentally at rather similar conditions 
by shock-compression experiments (27I 28 1. Moreover, 



the motion of charged particles trapped in such mag- 
netic fields causes the generation of radio waves. In fact, 
planetary radio experiments aboard the Voyager II flyby 
mission detected a wide variety of of radio emissions for 
both planets (29| . 

If large enough, amorphous or crystalline carbon clus- 
ters precipitate and sink towards the planetary center as 
sediment via gravitational settling. The corresponding 
release of energy has been estimated to be a substantial 
fraction of the internal heat production, which would ex- 
plain for instance why Neptune radiates more than twice 
the energy it receives from the sun Furthermore, it 
is also likely the cause for the externally observed high 
luminosity and could even contribute to the convective 
motions of its fluid interior of Neptune. The reason, why 
for Uranus no such internal heat flow mechanism could be 
established is still unknown Nevertheless, the similar- 
ity of the internal structures of these two planets suggests 
that the suppressed convection of Uranus may be a con- 
sequence of its closer proximity to the sun. In contrast, 
saturated hydrocarbons such as C2H6 and H2, being the 
products of the above ascertained decomposition of CH4 
at T = 4000 K, do not precipitate and instead rather rise 
to join the atmosphere. As a consequence, this process 
could be responsible for the anomalous abundance of H2 
in the atmospheres of both planets, and in the case of 
Neptune may also account for the observed wealth of at- 
mospheric C2H6, where it might be brought up from the 
deep interior by the just elucidated convection process. 
Therefore, the present results imply that deep chemical 
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processes such as phase transformations at extreme tem- 
peratures and pressures must be considered to model the 
interiors of giant gas planets more realistically. 

Even though, our calculations provide a consistent pic- 
ture of the deep chemistry of Neptune and Uranus, the 
remaining question is why no diamond formation could 
be observed, whereas experimentally it is reported to oc- 
cur from P = 20 GPa and T = 2000 K on. Due to the fact 
that liquid methane is optically transparent and can not 
simply be heated by a laser beam, it is therefore com- 
mon practice to include a noble metal absorber within 
laser-heated diamond anvil cell experiments. Spanu et al. 
reported that without a metallic absorber no formation 
of complex hydrocarbons and H2 at T = 2000 K could 
be determined, which not only agrees with the findings 
of the present work but also indicates that liquid CH4 
resides in a metastable state. On the contrary, at the 
pre sence of a nobel metal, liquid CH4 readily dissociates 

We conclude by noting that another possibility to ex- 
plain the discrepancy between theory and experiment 
may be the existence of a homogenous nucleation mech- 
anism, similar to the one recently proposed by Khali- 
ullin et al. for the direct graphite-to-diamond transition 

13, EH. 
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